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Abstract 

Laser-induced  breakdown  spectroscopy  (LIBS)  is  a  spectroscopic  technique 
where  output  from  a  pulsed  laser  is  focused  onto  a  target  in  order  to  create  an 
intense  plasma.  The  optical  emission  is  characteristic  of  the  elements  in  the  focal 
volume  and  can  be  used  for  elemental  analysis.  Research  on  the  detection  of 
nickel  in  solution  in  addition  to  solvent  detection  of  CCI4,  CHCI3,  C2CI4,  and 
C2HCI3  has  been  performed.  Breakdown  was  formed  at  the  sample  surface  via  a 
Q-switched  Nd:YAG  laser.  Initially,  operation  of  the  laser  was  at  1064 
nm/repetition  rate  of  5  Hz.  Experiments  were  also  performed  using  the  third 
harmonic  (355  nm)/repetition  rate  of  20  Hz.  Pulse  energy  was  maintained  at 
60  mJ.  The  spark  light  was  spectrally  resolved  and  detected  by  a  time-gated 
photodiode  array.  A  50  ps  gate  width/8  ps  time  delay  gave  detection  limits  of  56. 1 
mg/1  for  nickel  in  solution.  In  the  UV,  a  10  ps  gate  width/3  ps  time  delay  lowered 
detection  limits  down  to  29.4  mg/1.  Concentrations  spanned  from  50  to  1000  mg/1. 
Using  UV  excitation  (10  ps  gate  width/1  ps  time  delay),  saturated  solvent  solutions 
as  high  as  7.71  g/1  were  not  detectable. 
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Laser-Induced  Breakdown  Spectroscopy 
on  Solution  Samples  Using  Surface  Excitation 

I.  Introduction 

1.1  Overview 

Interest  in  monitoring  environmental  contamination  in  and  around  waste 
disposal  sites  and  research  facilities  has  grown  recently.  There  seems  to  be  an 
ever-increasing  need  to  quickly  and  accurately  measure  toxic-trace  elements  in  our 
environment.*  Several  methods  to  do  just  that  are  based  on  optical  spectroscopic 
techniques.  Such  examples  are  inductively-coupled  plasma  optical  emission 
spectroscopy  (ICP-OES),  graphite  furnace  atomic  absorption  spectroscopy 
(GF-AAS),  and  spark  discharge  optical  emission  spectroscopy  (SD-OES). 
Unfortunately,  these  methods  are  performed  in  a  laboratory  far  removed  from  the 
contaminated  site.  Samples  have  to  be  prepared  and  then  sent  out  to  a  laboratory 
for  analysis.  One  has  to  contend  with  a  time  delay  between  preparing  the  sample 
and  waiting  for  the  analytical  result  to  arrive  back  from  the  laboratory.^  With  the 
increasing  cost  of  sample  preparation,  the  demand  for  a  better  method  is  high. 

New  methods  are  needed  to  perform  these  concentration  measurements  in  situ 
while  still  maintaining  comparable  analysis  capability.  Laser-Induced  Breakdown 
Spectroscopy  (LIBS)  is  one  potential  method.  LIBS  uses  a  laser  to  vaporize. 


atomize,  and  excite  small  quantities  of  a  sample.  The  resulting  optical  emission 
can  then  be  used  for  elemental  analysis.* 

1.2  Problem  Foundation 

LIBS  is  of  interest  to  the  United  States  Air  Force  for  use  in  the 
environmental  monitoring  of  heavy  metals  (suspected  carcinogens)  in 
contaminated  groundwater.  With  this  in  mind,  this  study  will  begin  by  looking  at 
nickel  in  solution.  Additionally,  it  will  be  determined  whether  or  not  LIBS  is 
suitable  for  solvent  detection.  Specifically,  this  research  will  look  at  carbon 
tetrachloride  (CCI4),  chloroform  (CHCI3),  tetrachloroethylene  (C2CI4),  and 
trichloroethylene  (C2HCI3).  These  solvents  are  currently  being  used  in  Air  Force 
centers  worldwide  where  aircraft  are  cleaned  and  repainted.^  Unfortunately,  their 
toxicity  is  high  and  they  maybe  carcinogenic.'*  Inadequate  waste  disposal  of  these 
carcinogens  has  led  to  contaminated  groundwater  in  and  around  many  areas  of  the 
country.^  As  the  Department  of  Defense  (DoD)  tries  to  effectively  manage  safe 
handling  of  and  waste  disposal  of  these  carcinogens,  LIBS  may  potentially  allow 
for  easy  monitoring  of  these  activities.  Even  now,  under  a  Tri-Service  Reliance 
Agreement,  a  joint  U.S.  Army,  Navy,  and  Air  Force  effort  is  underway  to  develop 
and  test  optical  sensing  techniques  for  the  detection  of  soil  and  groundwater 
contaminants  of  interest  to  DoD. 
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1.3  Research  Objectives 

This  study  has  three  primary  objectives.  The  first  objective  is  to  construct 
an  apparatus  that  is  able  to  perform  Laser-Induced  Breakdovvn  Spectroscopy  on 
solution  samples.  The  second  objective  of  this  thesis  is  to  determine  if  LIBS  is 
able  to  detect  the  aforementioned  suspected  carcinogens  and  if  so;  the  final 
objective  of  this  study  is  to  generate  calibration  curves  of  analyte  signal  versus 
concentration  and  then  determine  detection  limits  for  these  samples. 

1.4  Sequence  of  Presentation 

Chapter  2  begins  with  some  of  the  background  leading  up  to  the  idea  of 
LIBS.  It  also  provides  some  of  the  theory  behind  LIBS.  Chapter  3  focuses  on  the 
experimental  apparatus  used  in  this  study.  It  also  discusses  sample  preparation 
techniques  used  in  this  research.  Chapter  4  presents  the  results  of  this  experiment 
while  Chapter  5  provides  a  discussion  of  the  data  inherent  to  this  study.  The  final 
chapter  reviews  the  objectives  of  this  thesis  and  offers  some  recommendations  for 
general  improvement  to  the  experimental  apparatus  as  well  as  some  conclusions 
representative  of  the  results. 


3 


II.  Literature  Review 


2.1  Background 

As  previously  mentioned,  ICP-OES,  GF-AAS,  and  SD-OES  are  currently 
being  utilized  for  elemental  analysis.  In  ICP-OES,  the  sample  is  prepared  as  an 
aerosol  which  is  then  injected  with  an  argon  gas  flow  into  a  heated  plasma  torch. 
Temperatures  of  the  resultant  plasma  are  between  6000  and  8000  K,  effectively 
atomizing  and  exciting  the  aerosol.  Multi-elemental  analysis  is  then  performed 
with  a  polychromator.  Detection  limits  are  usually  between  0.5  and  50  pg/l. 

On  the  other  hand,  GF-AAS  is  based  on  absorption  features  of  the  elements. 
The  sample  is  dissolved  and  a  droplet  of  the  solution  is  fed  into  a  graphite  furnace. 
The  furnace  effectively  atomizes  the  droplet  and  a  cloud  of  atoms  is  formed.  The 
absorption  of  radiation  of  an  element-specific  lamp  is  observed  to  determine 
element  concentration.  Unfortunately,  this  method  is  only  suitable  for  trace 
analysis  of  single  elements.  Limits  of  detection  for  this  technique  are  between 
0.001  and  30  pg/1. 

In  SD-OES,  a  ground  sample  is  atomized  and  excited  by  a  spark  discharge 
from  an  electrode.  Light  fi'om  the  resultant  plasma  is  observed  by  a  polychromator 
for  multi-elemental  analysis.  SD-OES,  however,  is  limited  to  only  substances  with 
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sufficient  electrical  conductivity  to  act  as  a  counter  electrode.  One  can  expect 
detection  limits  to  be  between  9  and  170  ppm.^ 

Notwithstanding  these  conventional  methods,  recent  interest  has  turned 
towards  the  laser  for  vaporization  and  excitation  of  the  target  sample.  The 
feasibility  of  using  lasers  as  excitation  sources  in  atomic  emission  spectroscopy 
dates  back  to  1962  when  Brech  and  Cross  first  demonstrated  this  capability.^ 

Since  then,  LIBS  has  been  used  in  the  analyses  of  gases,  liquids,  solids,  liquid  and 
solid  aerosols,  and  soil  samples.’  One  of  the  major  advantages  to  using  the  laser  is 
speed  and  simplicity.  Because  of  the  small  size  of  the  plasma  (around  20  mm  )  as 
well  as  its  low  energy  content  (<  500  mJ),  only  small  amounts  of  sample  are 
needed  for  elemental  analysis.  Depending  on  the  elements  present,  the  atoms 
excited  by  the  laser  may  be  neutral  or  ionized.  Species  are  identified  by  spectrally 
resolving  the  spark  light.  Other  advantages  over  conventional  methods  are; 

1)  No  auxiliary  analysis  equipment  is  needed  as  the  laser  vaporizes  and 
excites  the  sample  in  one  step. 

2)  The  absence  of  auxiliary  analysis  equipment  makes  this  method  very 
economical. 

3)  Since  LIBS  is  an  emission  technique,  direct  analysis  provides 
multielement  analysis  without  increased  complexity. 

4)  LIBS  is  useful  in  applications  requiring  noninvasive  analysis  because  the 
laser  spark  can  be  generated  in  remote  locations.* 
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Unfortunately,  there  are  a  couple  of  factors  that  can  make  the  use  of  LIBS 
difficult  as  an  analytical  technique.  One  drawback  is  the  possibility  of  inducing 
chemical  reaction  in  the  sample  under  analysis.  Laser  radiation  destroys  the 
molecular  bonds  of  the  pollutants  and  information  on  the  original  concentration 
should  be  determined  stoichometrically  by  studying  their  ionic  and  atomic  lines. 
The  atoms  of  the  molecules  tend  to  recombine  at  the  end  of  the  breakdown 
process,  producing  several  different  compounds  whose  concentrations  may  not  be 
easily  related  to  the  initial  concentration  of  the  pollutant.^  Secondly,  spectral  line 
intensities  last  for  several  microseconds  after  the  laser  pulse  and  then  decay  away 
in  a  complex  manner.  The  optimal  timing  for  spectral  acquisition  involves  a 
tradeoff  between  large  intensities  immediately  after  the  laser  pulse  and 
interferences  due  to  Stark  broadening  and  plasma  blackbody  radiation.^'* 

2.2  Theory  Behind  the  Experiment 

When  a  sufficient  amount  of  energy  from  an  energetic  laser  pulse  is  focused 
in  a  gas,  a  spark  plasma  is  formed  at  the  focus.  This  plasma  appears  as  a  bright 
flash  of  white  light  accompanied  by  a  loud  sound.  The  spectrally  dispersed  light 
emanating  from  the  spark  plasma  usually  exhibits  an  atomic  line  emission 
superimposed  on  a  broadband  continuum.  Typical  thresholds  for  gas  breakdown 
are  around  10^  W/cm^. 
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There  are  two  steps  leading  to  the  breakdown  of  this  gas.  The  first  step 
involves  the  generation  of  a  few  free  electrons  by  multiphoton  ionization  of  atoms 
and  molecules.  The  second  step  is  an  avalanche  ionization  of  matter  in  the  focal 
volume  to  form  a  plasma.  In  the  classical  sense,  electric  fields  of  the  optical  pulse 
accelerate  the  free  electrons  causing  them  to  collide  with  neutral  atoms.  This 
produces  an  isotropic  electron  energy  distribution.  Eventually,  the  electron  energy 
becomes  sufficiently  high  enough  to  collisionally  ionize  an  atom.  This  produces 
other  electrons  that  follow  the  same  course  and  cause  further  ionization.  This 
process  of  electron  multiplication  continues  throughout  the  pulse  and  results  in 
significant  breakdown  of  the  gas.  Following  this  breakdown,  the  plasma  expands 
radially  outward  in  all  directions  from  the  focal  volume.* 

The  laser  spark  can  also  be  generated  by  vaporizing  solid  surfaces. 

Focusing  a  high-powered  laser  onto  a  solid  sample  causes  a  complex  series  of 
events  to  occur.  The  intense  local  heating  experienced  by  the  sample  causes  the 
surface  temperature  of  the  target  to  rise  rapidly  at  a  rate  of  1 0’®  K/s.  As  heat  is 
conducted  to  the  interior  of  the  target,  a  thin  molten  layer  forms  just  below  the 
surface.  As  thermal  energy  deposited  at  the  surface  increases,  eventually  a  point  is 
reached  where  this  energy  exceeds  the  latent  heat  of  vaporization.  Once  this 
occurs,  heat  cannot  be  conducted  away  fast  enough  to  prevent  the  surface  from 
reaching  its  boiling  point  and  evaporation  occurs  at  the  surface. 
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Once  vaporization  occurs,  this  plasma  is  governed  by  gas  dynamic 
processes.  For  irradiances  just  above  threshold  (10^°  -  lO”  wW),  evaporation  at 
the  surface  proceeds  at  the  normal  boiling  point  of  the  material.  Conversely,  at 
irradiances  well  above  threshold  (10  W/m  ),  extremely  dense  plasmas  are 
formed.  Plasmas  produced  by  irradiances  between  the  two  extremes  are  usually 
the  most  conducive  to  LIBS  work.^* 

For  irradiances  within  the  specified  range  above,  vaporization  occurs 
rapidly  and  only  a  small  fi'action  of  the  sample  is  molten.  At  the  proper  wavelength 
(the  plasma  frequency),  the  vaporized  material  forms  a  luminous  plume  which 
extends  up  fi’om  the  surface  at  which  time  can  absorb  a  large  fi’action  of  the 
remaining  incident  pulse  energy.  Hence,  the  plume  effectively  shields  the  surface 
of  the  target  from  further  interaction  with  the  laser  pulse.  ^  After  removal  of  the 
laser  pulse,  the  plasma  continues  to  expand  radially  outward,  becoming  weaker 
and  diffuse  in  appearance. 

Hot  plasmas  are  rich  sources  of  excited  atoms  and  ions.  The  relative 
intensities  of  spectral  emission  lines  will  ultimately  lead  to  a  determination  of  the 
composition  of  the  target.**  One  of  the  most  important  requirements  for 
quantitative  analysis  is  the  concept  of  local  thermodynamic  equilibrium  (LTE). 
Simple  stated,  it  is  the  assumption  that  all  temperatures  in  the  plasma,  i.e.  ion, 
electron,  and  neutral,  are  all  the  same.  Moreover,  the  condition  of  LTE  means  that 
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rather  than  having  to  wait  for  the  whole  system  to  come  to  full  equilibrium,  it  is 
sufficient  to  characterize  a  localized  section  of  the  plasma  for  a  given  point  in  time. 
Assuming  that  LTE  has  indeed  been  established,  the  electron  temperature  defines 
the  excitation  state  of  the  plasma  and  elemental  concentrations  may  be  determined 
fi'om  the  intensities  of  the  spectral  lines. 

For  LTE  to  occur,  collisions  must  dominate  over  the  other  energy  transfer 
processes  such  as  radiative  decay  and  recombination.  This  will  establish  a 
Boltzmann-like  distribution  among  the  bound  energy  levels.  In  a  partially  ionized 
vapor,  collisions  will  be  the  dominant  mechanism  of  energy  transfer  because  of  the 
electron’s  high  velocities  and  their  long-range  coulomb  interaction.  However, 
radiative  processes  also  establish  the  populations  of  resonance  levels  of  the  major 
constituents  of  the  plasma.  Nevertheless,  the  small  concentrations  of  trace 
elements  present  in  the  plasma  ensure  that  electron  collisions  dominate  over 
radiative  processes  and  LTE  is  established  between  resonance  levels. 

Experimental  observation  of  emissions  from  laser-induced  plasmas  show 
that  initially,  an  intense  continuum  is  emitted  close  to  the  surface  of  the  target 
sample.  This  corresponds  to  the  emission  of  blackbody  (Bremsstrahlung)  radiation 
from  the  dense  plasma  and  typically  lasts  out  to  two  or  three  microseconds. 

Shortly  thereafter,  as  the  plasma  cools  and  expands  away,  line  emissions  tend  to 
dominate  the  spectrum,  with  the  first  ionized  species  being  emitted  close  to  the 
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target.  Atomic  lines  soon  appear  being  located  in  higher  regions  of  the  plume. 
Much  later,  emissions  from  simple  molecules  become  apparent.^* 

Even  though  gases  and  solids  have  been  fully  theoretically  explained,  there 
is  very  little  information  in  the  literature  that  describes  the  fundamental  physical 
processes  of  laser-induced  breakdovvn  in  liquids.  These  processes  are  not  very 
well  understood.  From  experiment,  it  follows  that  the  optical  breakdown  of  liquids 
differs  from  that  of  gases  in  its  external  manifestations.  Whereas  the  breakdown  of 
gases  results  in  a  single  laser  spark  in  the  lens  focus  (irrespective  of  focal  length), 
in  the  case  of  liquids,  long  focal  lengths  lead  to  the  production  of  several  laser 
sparks  along  the  laser  axis.  This  is  known  as  “multiplex”  breakdown.  In  fact,  this 
has  been  observed  in  gases  when  strongly  absorbing  microparticles  with  high 
concentration  were  introduced  into  the  laser  beam.  Because  of  this,  it  was  once 
thought  that  the  optical  breakdown  of  a  liquid  was  initiated  by  strongly  absorbing 
particles  suspended  in  it. 

Recently,  a  new  theory  about  optical  breakdown  not  associated  with  the 
idea  absorbing  particles  has  surfaced.  The  mechanism  considered  here  is  the  same 
as  that  in  gases,  namely,  electron  avalanche  development.  Such  a  process  is 
possible  if  there  are  gas  bubbles  in  the  focal  region  of  the  laser  beam.  In  a 
superheated  liquid,  gas  bubbles  can  be  generated  spontaneously.  The  development 
of  an  electron  avalanche  will  ultimately  lead  to  a  rise  in  pressure  inside  these 
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bubbles.  Consequently,  the  bubbles  expand  and  merge  into  one  large  bubble.  This 
is  known  as  stimulated  optic  coalescence  (SOC)  and  is  the  basis  for  this  new 
theory.  This  so-called  SOC  effect  is  also  accompanied  by  the  heating  of  gas 
particles  in  the  bubble  as  well  as  the  eventual  evaporation  of  its  liquid  wall.  What 
happens  is  that  the  large  bubble  continues  to  expand  at  the  expense  of  its  stored 
energy.  The  final  size  of  the  bubble  corresponds  to  the  size  of  “cavitation  bubbles” 
rising  to  the  surface  of  the  liquid  from  the  region  of  liquid  optical  breakdown.*^ 

2.3  Previous  Work  in  the  Field 

To  the  author’s  knowledge,  one  of  the  first  applications  of  LIBS  was  in 
1983.  Cremers  and  Radziemski  looked  at  the  detection  of  gases  in  air,  particularly 
chlorine  and  fluorine.  The  excitation  source  used  in  this  study  was  a  Nd:  YAG 
laser  operating  at  1 .06  pm.  Emission  from  the  plasma  was  spectrally  resolved  by  a 
Czemy-Tumer  spectrometer  and  monitored  by  a  thermoelectric  cooled 
photomultiplier  tube.  Time  resolution  of  the  plasma  light  was  performed  by  a 
boxcar  averager.  Detection  limits  were  established  at  8  ppm  for  chlorine  and  38 
ppm  for  fluorine. 

Another  study  in  1991  continued  to  investigate  the  application  of  LIBS  to 
gases  in  air.  Casini  et.  al.  looked  at  the  emission  spectra  of  a  variety  of  pollutants 
in  the  early  stages  of  plasma  development  (<  400  ns).  Again,  a  Nd:YAG  laser 
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operating  at  1 .06  pm  was  the  excitation  source.  Light  emitted  by  the  plasma  was 
focused  onto  the  slits  of  a  monochromator  and  the  resulting  signal  was  collected  by 
a  photomultiplier  tube.  Detection  limits  were  estimated  at  60  ppm  for  Cl  II,  200 
ppm  for  S  II,  200  ppm  for  P  H,  1 10  ppm  for  Na  H,  50  ppm  for  Hg  1, 1 30  ppm  for  Be 
II,  and  130  ppm  for  As  II. 

In  1994,  Lazzari  et.  al.  used  LIBS  to  measure  small  traces  of  heavy  metals  in 
the  atmosphere,  particularly  mercury.  The  experimental  setup  used  a  Nd:  glass 
laser,  pulse  width  8  ns,  to  deliver  intense  400  mJ  pulses.  Repetition  rate  of  the 
laser  was  at  10  Hz.  Laser  radiation  was  focused  onto  a  sample  contained  in  a  glass 
vessel  through  the  use  of  a  spherical  lens.  Light  from  the  plasma  was  sent  through 
an  optical  fiber  to  a  spectrometer  at  which  time  the  signal  from  the  spectrometer 
was  recorded  by  an  optical  multichannel  analyzer  (OMA).  To  attain  well-resolved 
spectra,  a  typical  averaging  of 400  samples  was  used.  To  obtain  precise 
determination  of  the  concentration  of  the  mercury  involved,  Lazzari  et.  al.  started 
out  with  a  small  amount  of  liquid  mercury  inside  an  experimental  cell.  After  a  few 
minutes,  an  equilibrium  concentration  of  gaseous  mercury  was  realized.  The 
liquid  mercury  was  then  removed  and  the  sample  was  ready  for  analysis.  The 
experiment  was  run  and  detection  limits  were  obtained.  Lazarri  et.  al.  found  that 
they  were  able  to  detect  concentrations  of  mercury  as  low  as  5  ppb.  This  is  quite 
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good  considering  the  fact  that  health  danger  thresholds  for  this  pollutant  is 
12  ppb.^ 

Finally  in  1996,  Robert  J.  Nordstrom  conducted  a  study  to  evaluate  the 
LIBS  spectral  characteristics  of  the  interference  from  N2  and  O2  gas  in  air.  A  CO2 
laser  was  used  as  the  excitation  source.  The  spectral  range  investigated  was  from 
350  to  950  nm.  Spectra  was  collected  by  a  grating  spectrometer  and  recorded  by  a 
photodiode  array.  Analysis  of  the  spectra  indicated  that  over  a  large  range  of  CO2 
laser  fluence,  all  spectra  could  be  coirelated  by  a  multiplication  factor.  This  meant 
that  changes  in  the  relative  spectral  contributions  from  N2  and  O2  were  not 
occurring.  Additionally,  it  indicated  that  the  contribution  of  different  atomic 
ionization  levels  in  the  spectra  did  not  change  appreciably.  This  result  is  important 
for  future  LIBS  work  where  air  interferences  must  be  considered.’^ 

In  addition  to  trace  elements  in  air,  LIBS  can  be  used  in  the  analysis  of 
solids.  In  1 985,  Belliveau  et.  al.  used  LIBS  to  detect  Cr  and  Mn  in  steel.  This 
study  employed  a  10  Hz,  1 .06  pm,  pulsed  Nd:YAG  laser  in  addition  to  two  echelle 
spectrometers  to  monitor  analyte  emissions.  Data  was  taken  and  calibration  curves 
were  developed.  Additionally,  a  correlation  study  between  dc  plasma  emission 
spectroscopy  and  LIBS  performed  on  Ni  showed  that  the  two  methods  gave  very 
similar  results.’^ 
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In  1 987,  Cremers  expanded  the  range  of  elements  to  which  LIBS  had  been 

applied  to.  He  investigated  Mo,  Ni,  Al,  Ti,  Cu,  Fe,  Pb,  and  In.  The  source  of 

excitation  for  this  experiment  was  a  Q-switched  Nd:YAG  laser.  This  time, 

however,  light  was  collected  by  means  of  a  fiber  optic  cable  and  transmitted  to  a 

spectrograph  for  spectral  resolution.  Spectra  was  recorded  by  a  diode  array  rapid 

scanning  spectrometer  (DARSS).  This  work  demonstrated  that  metals  could  be 

detected  via  LIBS  at  up  to  distances  of  2.4  m  between  the  sample  and  the  focusing 

1 8 

lens-light  collection  optics. 

Work  by  Grant  et.  al.  proved  to  be  similar  in  nature  to  Cremers.  In  1 991, 
they  used  LIBS  to  perform  elemental  analysis  on  iron  ore.  Unlike  previous 
experiments,  an  excimer  laser  was  the  excitation  source.  Again,  a  fiber  optic  cable 
was  used  to  collect  and  transmit  the  light  emitted  by  the  plasma  to  a 
monochromator.  A  photomultiplier  tube  recorded  the  spectra.  Calibration  curves 
were  generated  by  plotting  the  ratio  of  the  intensity  of  each  element  (Ca,  Si,  Mg, 
Al,  and  Ti)  to  that  of  iron  ore.  Detection  limits  were  on  the  order  of  0.01%.^ 

In  1 994,  Alexander  et.  al.  set  out  to  measure  detection  limits  for  chromium. 
The  experiment  used  a  Q-switched  Nd.  YAG  laser  operating  at  1064  nm  with  a 
pulse  width  of  10  ns  as  an  excitation  source.  Light  Ifom  the  laser  was  focused  onto 
the  sample  by  a  plano-convex  BK7  lens.  Next,  light  emitted  by  the  plasma  was 
focused  into  an  optical  fiber  for  transmission  to  an  OMA  consisting  of  a 
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spectrometer  and  a  photodiode  array.  They  were  able  to  generate  a  calibration 
curve  for  chromium  with  a  detection  limit  down  to  100  ppb.^ 

That  same  year,  Thiem  et.  al.  used  LIBS  in  a  ultra-high  vacuum  to  perform 
elemental  analysis  on  NIST  transition  metal  alloys.  The  second  harmonic 
(k  =  532  nm)  of  a  Nd;YAG  laser  was  used  as  the  excitation  source  to  ablate  the 
surface  of  a  standard  metal  sample.  The  sample  was  kept  on  an  X-Y  translation 
table  that  was  constantly  moving  throughout  the  course  of  irradiation.  Light 
emitted  by  the  plasma  was  collected  by  means  of  a  fiber  optic  bundle  which  then 
guided  the  light  to  the  spectrograph  and  photodiode  array.  Data  was  taken  and 
detection  limits  were  established:  0.08%  for  Al,  1 .45  %  for  Cr,  0.70%  for  Co, 
0.01%  for  Cu,  0.54%  for  Fe,  0.07%  for  Mn,  0.001%  for  Ni,  0.04%  for  Si,  0. 1 1% 
for  Ti,  and  0.44%  for  Zn.^^ 

Two  years  later,  Pakhomov  et.  al.  used  LIBS  to  detect  lead  in  concrete.  The 
breakdown  was  formed  at  the  sample  surface  by  a  pulsed  Nd:YAG  laser  operating 
at  1 .06  pm.  Light  fi’om  the  plasma  was  imaged  onto  a  quartz  optical  fiber.  The 
opposite  end  of  the  fiber  was  attached  to  the  entrance  slit  of  a  Czemy-Tumer 
spectrograph.  Light  was  detected  by  a  gated  photodiode  array.  Contamination 
levels  were  inferred  Ifom  the  ratio  of  the  emission  line  of  lead  to  a  known 
reference  line.  Lead  contamination  was  determined  down  to  10  ppm.^° 
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Recent  interest  has  turned  towards  the  analysis  of  liquids  via  LIBS. 

Cremers  et.  al.  first  looked  at  LIBS  as  a  way  to  analyze  liquids  in  1984.  They 
looked  at  the  dielectric  breakdown  of  several  elements  in  water.  They  noted  that 
this  breakdown  occurs  with  focused  laser  powers  of  10^®  - 10^^  W/cm^  as 
compared  to  that  of  10^-10^  W/cm^  in  air.^’  However,  this  threshold  value  is 
affected  by  such  factors  as  impurities  contained  in  the  medium,  irradiation 
geometry,  and  laser  pulse  width.  Output  fi'om  a  Nd:YAG  laser  was  focused  into 
various  samples  by  a  pair  of  5  cm  focal  length  lenses.  The  pulses  were  focused 
into  the  cell  fi'om  the  side  because  pulses  entering  from  the  top  intercepted  bubbles 
from  previous  sparks  thereby  interfering  with  spark  formation.  Light  emitted  by 
the  plasma  was  spectrally  resolved  and  then  detected  by  a  photomultiplier  tube 
(PMT).  Detection  limits  were  established  for  Li,  Na,  K,  Rb,  Cs,  Be"^,  Mg"^,  Ca"^, 
and  B.  For  example,  limits  of  detection  for  the  alkali  metals  ranged  from  0.006  to 
1 .2  pg/ml.  On  the  other  hand,  detection  limits  for  the  nonalkali  elements  proved  to 
be  much  higher,  anywhere  from  1  to  1200  |ig/ml.^* 

In  1987,  Wachter  et.  al.  looked  at  uranium  in  solution.  They  note  that 
previous  work  analyzed  liquids  by  the  formation  of  the  laser  spark  in  a  bulk  liquid. 
This  avoided  the  problem  of  spark  perturbation  due  to  surface  agitation  and 
splashing  of  the  liquid,  but  resulted  in  a  lower  temperature  and  hence  reduced 
excitation  capability  as  that  compared  to  a  gas.  Previous  work  has  shown  that  with 
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bulk  excitation,  uranium  concentrations  as  high  as  300  g/1  could  not  be  detected. 
However,  with  surface  excitation,  uranium  was  detected  at  concentrations  at 
several  orders  of  magnitude  below  300  g/1.  Detection  limits  are  improved  with 
surface  excitation  because  the  solution  is  vaporized  into  hot  plasma  formed  above 
the  liquid  surface  where  the  sample  is  more  efficiently  excited.^^  Other  advantages 
associated  with  surface  excitation  are  that  light  emitted  by  the  plasma  will  not  be 
absorbed  by  other  substances  in  the  solution  and  that  particles  suspended  in  the  air 
will  not  perturb  spark  formation  nor  scatter  plasma  light  to  any  great  extent.^"* 
Wachter  et.  al.  focused  pulses  from  a  Q-switched  Nd:YAG  laser  onto  the  surface 
of  a  solution  thereby  generating  the  laser  spark.  The  plasma  light,  viewed  through 
the  side  of  the  glass  vial  containing  the  solution,  was  then  focused  onto  the  slit  of  a 
spectrograph.  The  spectrally  resolved  light  was  then  recorded  by  a  DARSS 
detector.  Using  surface  excitation,  Wachter  et.  al.  were  able  to  lower  detection 
limits  for  uranium  down  to  0. 1  g/1. 

In  1991,  Vlasov  et.  al.  studied  the  question  of  how  the  principal  elements  of 
salt  water  solutions  change  in  time  as  a  result  of  breakdown  at  the  water  surface. 
Their  principal  focus  was  on  estimating  the  detection  limits  and  accuracy  to  which 
trace  elements  present  in  sea  water  could  be  determined.  A  Nd:YAG  laser  was  the 
excitation  source.  Spectra  obtained  were  either  time  integrated  or  taken 
instantaneously.  Comparison  of  the  two  show  that  with  time,  the  continuum 
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background  diminishes  rapidly  while  the  contrast  of  the  line  spectrum  increases. 
For  example,  the  N 11  lines  emit  only  during  the  first  few  hundred  nanoseconds. 
This  means  that  for  all  practical  purposes,  the  atmospere  does  not  penetrate  into 
the  plume  and  only  the  surface  layer  of  air  contributes  to  the  emission  from  the 
plume.  They  also  note  that  the  line  reaches  maximum  intensity  within  2  ps  and 
virtually  disappears  after  6  ps.  However,  Na  remains  strong  even  out  to  8  ps.  A 
calibration  curve  for  Na  reveals  a  linear  dependence  for  concentrations  between 
0.004  -  0.2  %  by  mass.^^ 

A  few  years  later,  Cespedes  et.  al.  looked  at  the  analysis  of  liquids  via  LIBS. 
A  Nd.  YAG  laser  =  1064  nm)  operating  at  10  Hz  with  a  pulse  duration  of  9  ns 

was  used  as  the  excitation  source.  Light  from  the  plasma  was  collected  by  a  lens 
and  imaged  onto  the  entrance  slit  of  a  0.64  m  spectrometer.  The  dispersed  light 
exited  the  spectrometer  and  was  recorded  by  an  intensified  photodiode  array. 
Aqueous  solutions  of  Cd,  Co,  Cr,  Cu,  Fe,  Mn,  Ni,  Pb,  Tl,  Zn,  Ag,  and  Hg  were 
prepared.  Data  was  taken  and  detection  limits  were  obtained;  70  ppm  for  Mn,  300 
ppm  for  Pb,  6  ppm  for  Tl,  100  ppm  for  Cr,  and  1500  ppm  for  Fe.  All  other 
solutions  were  not  detectable  even  at  very  high  concentrations  (10,000  ppm).^^ 

In  1 996,  Knopp  et.  al  investigated  LIBS  as  an  analytical  tool  for  the 
detection  of  metal  ions  in  solution.  An  excimer  (k  =  308  nm)  pumped  dye-laser 
system  was  the  excitation  source.  The  dye-laser  pulse  was  at  500  nm  with  an 
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output  of  22  ±  2  mJ/pulse.  Plasma  emitted  light  was  monitored  by  a  polychromator 
and  detected  by  an  intensified  time  gated  photodiode  array.  Detection  limits  for 
the  metal  ions  were  500  mg/1  for  Cd^^,  12.5  mg/1  for  Pb^"^,  6.8  mg/1  for  Ba^^,  0. 13 
pg/1  for  Ca^"^,  13  pg/1  for  Li"^,  and  7.5  pg/1  for  Na"^.  Hg^"^  and  Er^"*^  were  not 
detectable  in  solution. 

To  the  author’s  knowledge,  there  has  not  been  any  research  on  solvent 
detection  via  LIBS.  Solvents  pose  a  great  health  hazard  to  humans  and  this  area  of 
work  needs  to  be  examined.  This  study  will  answer  the  question  of  whether  or  not 
LIBS  is  a  viable  method  for  solvent  detection. 
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III.  Experimental  Design  and  Methodology 


3.1  Instrumentation 

The  experimental  setup  used  in  this  study  is  shown  in  Figure  3.1.  The 
output  from  a  Nd:YAG  laser  (Spectra-Physics  Quanta  Ray  DCR-3)  operated  in 
Q-switched  mode  was  focused  via  a  10  cm  plano-convex  lens  onto  the  surface  of  a 
liquid  sample.  A  prism  was  used  to  direct  the  laser  pulse  downward  into  a  cuvette 
where  the  solution  sample  was  held.  Initially,  the  laser  was  run  at  1064  nm  with  a 
repetition  rate  of  5  Hz.  But  soon  the  choice  was  made  to  run  the  laser  at  the  third 
harmonic  (A,  =  355  nm).  Repetition  rate  was  now  at  20  Hz  (unless  otherwise 
noted).  The  purpose  of  doing  this  was  to  look  at  wavelength  dependencies  on 
LIBS  signal  as  well  as  detection  limits.  In  both  cases,  the  laser  was  set  to  deliver 
60  mJ/pulse  (unless  otherwise  stated).  This  meant  that  the  laser  fluence  on  target 
was  approximately  10^  - 1  O'*  J/cm^.^  Not  shown  in  Figure  3.1  is  a  pyroelectric 
detector  (Scientech,  Model  380101).  This  was  used  to  monitor  pulse  energy. 

Pulse  width  of  the  laser  was  about  15  ns.  An  average  of  75  laser  pulses  (50  at 
1064  nm)  was  used  to  obtain  sufficient  signal-to-noise  (S/N)  ratios. 

Fine  tuning  of  the  laser  focus  onto  the  surface  of  the  solution  sample  was 
accomplished  through  the  use  of  2-dimensional  translation  stages.  Both  the  lenses 
(plano-convex  and  convex)  as  well  as  the  cuvette  holder  were  mounted  onto  these 
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Figure  3.1  Schematic  of  Experimental  Apparatus 

translation  stages.  After  completion  of  the  adjustment,  the  slits  of  the  spectrograph 
(Acton  Research  Corporation,  Spectra-Pro  275)  were  aligned  with  respect  to  the 
location  on  the  liquid  surface  where  the  laser  was  focused. 

The  spectrograph  was  oriented  orthogonally  to  the  laser  axis  in  order  to 
collect  light  emitted  fi'om  the  focal  volume.  A  slit  width  of  100  pm  was  used  for 


nickel  in  solution  while  50  pm  was  used  for  solvent  detection  (unless  otherwise 
noted).  Wavelength  region  dictated  what  grating  was  to  be  used.  The  1200 


grooves/mm  grating  was  blazed  at  300  nm  and  was  used  for  wavelengths  less  than 
400  nm.  The  600  grooves/mm  grating  was  blazed  at  500  nm  and  was  used  in  the 
wavelength  region  of 400-760  nm.  It  was  determined  that  the  spectral  resolution 
for  the  1200  grooves/mm  grating  was  0.6  nm  while  it  was  1  nm  for  the  600 
grooves/mm  grating.  Once  the  spectrograph  dispersed  the  light,  the  intensified 
photodiode  array  (EG&G  PARC,  Model  1420)  would  record  this  light  as  an  analog 
signal.  The  detector  interface  (EG&G  PARC,  Model  1461)  then  would  take  this 
signal  and  convert  it  into  a  digital  one  where  it  the  would  then  be  displayed  on  a 
computer  monitor  as  intensity  versus  pixel  number. 

A  wavelength  calibration  was  used  to  convert  the  pixel  number  scale  to  one 
of  a  wavelength  scale.  This  was  done  through  the  use  of  Ne,  Hg,  Ar,  Kr,  and  Xe 
penlamps.  It  involved  setting  the  spectrograph  to  a  predetermined  wavelength  and 
capturing  light  from  the  penlamps.  Spectra  again  was  displayed  as  intensity  versus 
pixel  number.  Wavelengths  of  these  elements  are  well  known  and  were  used  to 
match  the  peaks  in  the  spectra  with  known  wavelengths.^^  A  linear  fit  of 
wavelength  to  pixel  number  produced  a  good  calibration  for  the  spectra.  This  was 
done  every  20  nm  (for  both  gratings).  Calibration  error  was  up  to  0.2  nm  for  the 
1200  grooves/mm  grating  while  it  was  up  to  0.5  nm  for  the  600  grooves/mm 
grating.  Before  all  spectra  were  subject  to  analysis,  the  background  generated  by 
noise  in  the  photodiode  array  and  associated  electronics  would  be  recorded  and 
then  subtracted  from  the  spectra. 
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So  far,  the  FND-100  detector  (EG&G  PARC),  the  oscilloscope  (LeCroy, 
Model  1420),  the  digital  delay  (Stanford  Research  Systems,  Model  DG535),  and 
the  pulse  amplifier  (EG&G  PARC,  Model  1304)  have  not  been  discussed.  These 
were  components  of  the  equipment  used  in  the  timing  of  the  system.  The  FND- 
100  detector  was  used  to  detect  the  laser  pulse  and  initiate  timing  of  the  system. 
The  oscilloscope  was  used  to  monitor  this  timing.  It  displayed  voltage  versus  time. 
The  digital  delay  generated  a  TTL  (transistor-transistor  logic)  pulse  that  triggered 
the  detector  interface  and  associated  electronics.  Finally,  the  pulse  amplifier  sent  a 
gating  pulse  to  the  photodiode  array.  Timing  of  the  system  is  shown  in  Figure  3.2. 
This  is  a  trace  from  the  oscilloscope.  As  can  be  seen,  the  laser  started  at  time  t  =  0. 
Ninety  nanoseconds  later,  the  digital  generator  generated  a  TTL  pulse.  One 
hundred  and  eighty  nanoseconds  after  initiation  of  the  laser  trigger,  the  detector 
interface  received  that  TTL  pulse.  Two  hundred  and  twenty  nanoseconds  later,  the 
pulse  amplifier  was  triggered  at  which  time  the  photodiode  array  was  sent  the 
gating  pulse.  The  photodiode  array  then  took  its  measurements  at  which  time  the 
system  waited  for  the  cycle  to  begin  again.  When  the  laser  was  operated  at  5  Hz,  a 
cycle  was  completed  every  200  ms,  while  for  20  Hz,  a  cycle  took  50  ms. 

One  also  needed  to  choose  a  gate  width  and  a  time  delay  that  would  give 
rise  to  maximized  S/N  ratios.  To  do  this,  the  experiment  was  run  several  times 
using  a  solution  sample  (1  g^  for  Ni  and  0.971  g/1  for  CCI4)  at  different  gate 
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Time  (ns) 


Figure  3.2  Timing  of  the  System 


widths  and  time  delays.  After  looking  over  the  results,  a  50  fis  gate  width  and 
a  8  )rs  time  delay  was  chosen  for  use  at  the  fundamental  wavelength  (1 .06  |am). 
However,  for  the  third  harmonic  (k  =  355  nm),  a  10  |is  gate  width  and  a  3  fj.s  time 
delay  was  used.  For  solvent  detection  (CCI4),  a  1  (is  time  delay  was  utilized. 
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3.2  Sample  Preparation 

In  order  to  generate  a  calibration  curve,  one  must  run  the  experiment 
several  times  at  differing  concentrations.  To  accomplish  this  task,  nickel  in 
solution  was  diluted  according  to  the  following  formula; 


Molaritybefore  X  Volumebefore  ^  Molarity  after  X  Volume  after  (3.1) 


where  molarity  (M)  is  in  units  of  moles/liter. 

For  example,  if  one  had  a  solution  with  an  initial  concentration  of  12  M,  how  much 
water  must  be  added  to  decrease  the  molarity  of  the  new  solution  down  to  0. 15  M 
and  obtain  500  ml  of  the  new  batch?  Using  the  formula  above,  one  would  need  to 
measure  out  6.25  ml  of  the  12  M  solution  and  add  enough  water  to  the  solution 
until  the  total  volume  reached  500  ml.^^ 

Solvents,  on  the  other  hand,  were  prepared  from  concentrates.  First  one 
needed  to  determine  their  saturation  limit,  i.e.,  how  much  of  the  concentrate  would 
the  water  absorb.  The  following  formulas  were  used.  Conditions  were  at  25°C 
and  units  were  in  moles/liter  (M): 

For  CHCI3,  log  =  - 1.19  (3 .2) 

Cw  =  0.0646  M  =  7.71  g/1 

ForCCU,  log  2.20  (3.3) 

Cw  =  0.00631  M  =  0.971  g/1 
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ForC2HCl3, 


log  Cw  =  -  2.04 


(3.4) 


Cw  =  0.00912  M=  1.20  g/1 

ForC2Cl4,  log  -3.04  (3.5) 

Cw  =  0.151  g/1 

One  would  prepare  the  base  (saturated)  solutions  by  first  measuring  out  the  proper 
amount  of  concentrate  needed.  A  pan  balance  was  used  to  perform  this  task. 

Next,  a  liter  of  water  was  weighed  out  (density  ~  1000  g/1)  and  then  added  to  the 
solvent.  Once  the  base  solutions  were  available,  the  dilutions  were  prepared  in  the 
same  manner.^^ 
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rv.  Results 


4.1  Nickel  in  Solution  (IR  Excitation  at  1064  nm) 

A  typical  spectrum  covering  the  wavelength  region  from  335  to  385  nm  for 
nickel  in  solution  is  shown  below  in  Figure  4. 1 .  Spectral  line  assignments  are 
shown  in  Table  4. 1  (Transitions  shown  are  from  lower  state  to  upper  state.)  All 
nickel  lines  were  atomic  in  nature.  Ions  were  not  detected  in  this  spectrum.  As 
one  can  see,  there  are  many  discernible  peaks  to  choose  from  for  elemental 
analysis.  In  this  case,  the  352.45  nm  line,  labeled  “D”  in  the  figure,  was  chosen. 
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Figure  4.1  Ni  (1  g/1  -  IR)  at  50  ps  Gate  Width  and  8  ps  Time  Delay 
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Table  4.1  Spectrum  Line  Assignments  for  Figure  4.1 


Letter 

Wavelength  (nm) 

Transition 

A 

341.48  nm 

B 

345.85  nm 

C 

349.30  nm 

D 

352.45  nm 

E 

356.64  nm 

F 

361.94  nm 

my^ESESlil 

This  was  the  logical  choice  as  it  had  a  very  high  intensity  as  that  compared  to  the 
others.  Also  note  that  water  plays  an  integral  part  in  it’s  contribution  to  the 
spectrum.  This  can  be  seen  from  Figure  4.2.  However,  only  one  peak  is  clearly 
discernible  in  the  water  spectrum.  It  is  due  to  NH  (A  -  X  "  system).^^’^'* 

Note  that  as  one  decreases  the  time  delay,  the  continuum  background  due  to 
Bremsstrahlung  radiation  becomes  the  dominant  feature.  This  is  shown  in  Figures 
4.3  and  4.4.  At  the  beginning  of  the  experiment,  the  optimal  conditions  were  not 
known.  Using  a  pulse  amplifier,  spectra  were  gated  and  observation  times  were 
chosen.  This  was  done  to  maximize  S/N  ratios.  Since  the  optimum  time  is 
somewhat  dependent  on  the  species  involved,  it  should  be  determined  for  each 
species  investigated.  The  spectra  were  recorded  and  then  a  S/N  ratio  was 
determined  for  each  plot.  The  signal  can  be  determined  in  one  of  many  different 
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ways.  One  could  use  the  value  of  the  intensity  at  the  peak  height.  But  this  does 
not  take  into  account  how  much  of  the  intensity  is  due  to  the  continuum 
background  nor  does  it  consider  the  variation  of  the  background  as  a  function  of 
concentration.  It  is  much  better  to  use  the  difference  between  peak  and  base 
heights  as  a  measure  of  the  signal.  Alternatively,  one  could  use  the  area  under  the 
curve  as  a  measure  of  the  signal.  The  problem  with  this  method  is  where  does  the 
spectral  line  start  and  the  continuum  background  end.  One  must  be  very  careful  in 
choosing  which  spectral  line  to  use  for  elemental  analysis.^  The  noise  was 


c 

u 

ea 


a 

u 


800 


600 


400 


200 


0 


-200 


320 


340 


360 


380 


400 


Wavelength  (nm  ) 


Figure  4.2  Deionized  Water  (IR)  at  50  jis  Gate  Width  and  8  ps  Time  Delay 
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estimated  by  calculating  the  standard  deviation  of  the  continuum  background 
(RMS  noise)  in  a  portion  of  the  spectrum  that  was  essentially  featureless.  The 
value  of  the  intensities  in  that  region  were  known  so  that  the  standard  deviation 
was  easily  obtainable.  In  comparing  S/N  ratios  with  area/noise  ratios,  one  obtained 
a  much  higher  value  for  the  S/N  ratio.  This  is  shown  in  Appendix  A.  Hence,  it 
was  easily  justifiable  to  use  the  differencing  technique  as  the  measure  of  the  signal 
in  this  research.  It  can  be  seen  fi-om  Figure  4.5  that  the  trend  shows  that 
somewhere  between  6  and  10  ps  gives  maximized  S/N  ratios.  A  value  of  8  ps  was 
chosen  as  the  optimal  delay  time.  Note  that  the  curve  follows  a  quadratic  nature. 
The  statistical  error  is  also  plotted  as  error  bars  in  the  figure.  (Another  source  of 
error  in  the  experiment  is  systematic  error.  This  is  briefly  discussed  in  Appendix 
B.)  Differing  gate  widths  were  also  examined.  There  came  a  point  when  the  S/N 
ratio  dropped  as  one  increased  the  gate  width.  This  was  probably  due  to  increased 
noise  as  the  gate  width  was  increased.  Table  4.2  summarizes  the  different  time 
delays  and  gate  widths  that  were  looked  at.  The  calibration  curve  for  nickel  in 
solution  was  taken  by  monitoring  the  atomic  line  at  352.45  nm  8  ps  after  plasma 
initiation.  It  is  shown  in  Figure  4.6.  Diluted  solutions  spanned  from  100  to  1000 
mg/1.  The  curve  had  a  correlation  coefficient  of  0.999.  Detection  limits  were 
established  by  the  following  formula: 
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Cl=UL 

S 


(4.1) 


where  N,  as  before,  is  the  standard  deviation  of  the  continuum  background  and  S 
is  the  slope  of  the  calibration  curve.^  Slope  of  the  curve  was  0.556 1/mg  while  the 
standard  deviation  of  the  background  (from  the  100  mg/1  solution)  was  10.4. 
(Noise  levels  at  the  other  concentrations  are  shown  in  Appendix  C.)  Detection 
limits  for  nickel  in  solution  using  IR  excitation  were  about  56. 1  mg/1. 


Figure  4.5  S/N  Ratio  Versus  Time  Delay  at  352.45  nm 
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Table  4.2  Results  from  Gate  Width  Experiment  for  Ni  (IR) 


Gate  Width  (ps) 

Time  Delay  (ps) 

Signal 

Noise 

S/N  Ratio 

50 

0 

904 

118 

7.66 

50 

0.5 

569 

61.8 

9.21 

50 

1 

902 

52.8 

17.1 

50 

4 

409 

16.3 

i 

25.1 

50 

8 

561 

8.38 

66.9 

50 

12 

202 

6.72 

30.1 

1 

8 

40.0 

8.22 

4.87 

10 

8 

122 

9.58 

12.7 

25 

8 

332 

13.4 

24.8 

75 

8 

341 

9.11 

37.4 

4.2  Nickel  in  Solution  (UV  Excitation  at  355  nm) 

A  typical  spectrum  covering  the  wavelength  region  from  335  to  385  nm  for 
nickel  in  solution  is  shown  below  in  Figure  4.7.  Spectral  line  assignments  are 
shown  in  Table  4.3.  Again,  analytical  lines  were  atomic  in  nature.  No  ions  were 
detected.  Note  that  the  intensity  of  these  nickel  lines  are  much  greater  than  before 
-  by  a  factor  of  3  to  1 .  In  addition  to  the  352.45  nm  line,  analytical  lines  at  341 .48 
and  361 .94  nm  were  also  chosen  for  elemental  analysis.  This  was  done  because  of 
possible  interferences  from  the  355  nm  laser  line.  Once  more,  the  continuum 
background  due  to  Bremsstrahlung  radiation  became  the  dominant  feature  as  the 


Figure  4.7  Ni  (1  g/I  -  UV)  at  10  ps  Gate  Width  and  3  ps  Time  Delay 
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Table  4.3  Spectrum  Line  Assignments  for  Figure  4.7 


Letter 

Wavelength  (nm) 

Transition 

A 

341.48  nm 

IQjgBgllllll 

B 

345.85  nm 

C 

349.30  nm 

mmmamm 

D 

352.45  nm 

Kgiofoomi 

E 

361.94  nm 

time  delay  was  decreased.  This  can  be  seen  in  Figure  4.8.  Again,  spectra  were 
gated  in  order  to  determine  an  optimal  delay  time.  This  is  shown  in  Figure  4.9. 
The  trend  clearly  shows  that  somewhere  around  3  ps  gave  maximized  S/N  ratios. 
As  before,  the  curve  had  a  quadratic  nature.  Differing  gate  widths  were  also 
looked  at.  Table  4.4  summarizes  these  different  conditions.  (Data  taken  at  0.5  ps 
time  delay  was  at  10  Hz  and  28  mJ.)  Looking  at  Tables  4.2  and  4.4,  one  can  see 
that  there  is  no  real  wavelength  dependence  on  the  S/N  ratios.  Even  though  there 
was  an  increase  in  LIBS  signal  using  UV  excitation,  noise  had  increased  as  well. 
The  calibration  curve  for  nickel  in  solution  taken  by  monitoring  the  atomic  line  at 
361 .94  nm  3  ps  after  plasma  initiation  is  shown  in  Figure  4.10.  Diluted  solutions 
spanned  from  50  to  1000  mg/1.  The  curve  had  a  correlation  coefficient  of  0.994 
while  slope  of  the  curve  was  1 .52  1/mg.  The  RMS  noise  (from  the  50  mg/1 
solution)  was  20.0.  (Noise  levels  at  the  other  concentrations  are  listed  in 
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Intensity  (arb.  units) 


Figure  4.8  Ni  (1  g/1  -  UV)  at  10  jis  Gate  Width  and  1  ns  Time  Delay 


Figure  4.9  S/N  Ratio  Versus  Time  Delay  at  361.94  nm 
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Intensity  (arb.  units) 


Table  4.4  Results  from  Gate  Width  Experiment  for  Ni  (UV) 


Gate  Width  (^is) 

Time  Delay  (ps) 

Signal 

Noise 

S/N  Ratio 

10 

1 

5318 

186 

28.6 

10 

3 

1416 

16.1 

88.0 

10 

6 

367 

8.11 

45.3 

10 

9 

254 

16.1 

15.8 

10 

12 

95 

10.7 

8.88 

3 

0.5 

282 

21.2 

13.3 

7 

0.5 

324 

17.4 

18.6 

10 

0.5 

256 

10.6 

24.2 

14 

0.5 

244 

10.8 

22.6 

18 

0.5 

226 

11.0 

20.5 

Figure  4,10  Calibration  Curve  for  Ni  (UV)  at  361.94  nm 
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Appendix  C.)  The  detection  limit  for  this  curve  using  UV  excitation  was  about 
39.5  mg/1.  Data  for  all  three  calibration  curves  is  shown  in  Table  4.5.  The  average 
detection  limit  was  about  33.6  mg/1  with  a  standard  deviation  of  5.26  mg/1. 


Table  4.5  Calibration  Curve  Data  for  Ni  (UV) 


Wavelength 

Correlation  Coef. 

Slope  of  Curve 

Detection  Limit 

341.48  nm 

0.988 

1.88 

31.9  mg/1 

352.45  nm 

0.984 

2.04 

29.4  mg/1 

361.94  nm 

0.994 

1.52 

39.5  mg/1 

Average 

— 

— 

33.6  +  5.26  mg/1 

4.3  Pure  Solvent  Detection  (UV  Excitation) 

Figure  4.1 1  is  a  spectrum  for  CCI4  covering  the  wavelength  region  of  705 
to  805  nm.  Two  distinct  chlorine  lines  can  be  seen  at  725.66  and  754.71  nm.35,36 
They  are  atomic  in  nature  and  ions  were  not  detected.  Additionally,  second  order 
cyanogen  lines  (violet  system,  B  2^  -  X  ^S)  as  well  as  nitrogen  lines  were  also 
discovered.34, 37,38,39  xhg  nitrogen  lines  appear  from  the  diffusion  of  air  into  the 
solution  sample  from  the  surface  layer  of  air  above  it.  Line  assignments  for  this 
spectrum  are  listed  in  Table  4.6.  A  similar  spectrum  is  also  observed  for  CHCI3, 
C2CI4,  and  C2HCI3  (see  Appendix  D).  One  thing  that  was  noticed,  however,  was 
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Intensity  (  x  10^  arb.  units) 


(nm) 


805  nm 


re  4.11 


Letter 


Element  Wavelength  Transition 


CN  (2nd  order) 

359.04  nm 

Cl  I 

725.66  nm 

NI 

742.36  nm 

NI 

744.23  nm 

NI 

746.83  nm 

Cl  I 

754.71  nm 

CN  (2nd  order) 

385.47  nm 

CN  (2nd  order) 

386.19  nm 

CN  (2nd  order) 

387.14  nm 

CN  (2nd  order) 

388.34  nm 
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the  fact  that  the  intensities  for  the  chlorine  bands  were  not  all  the  same.  For 
example,  it  can  be  seen  from  Table  4.7  that  single  bond  molecules  have  higher  Cl 
line  intensities  than  their  double  bond  counterpart,  i.e.,  CCI4  >  C2CI4  and  CHCI3  > 
C2HCI3.  Additionally,  chlorocarbons  have  higher  Cl  line  intensities  than  their 
hydrocarbon  counterpart,  i.e.,  CCI4  >  CHCI3  or  C2CI4  >  C2HCI3. 

Figure  4. 12  shows  a  CCI4  spectrum  covering  the  wavelength  region  of  585 
to  685  nm.  Only  one  peak  tends  to  stand  out  due  to  the  large  continuum 
background.  Surprisingly,  it  is  a  hydrogen  peak  (at  656.29  nm).^^’^^  But  how 
could  this  be?  Checking  the  bottle  the  CCI4  was  stored  in,  it  was  noted  that  there 
were  10  ppm  of  water  in  the  solution.  Mystery  solved!  A  similar  peak  is  observed 
in  the  other  solvents  as  well  (see  Appendix  D). 

Figure  4. 1 3  is  a  C2CI4  spectrum  covering  the  wavelength  region  of 465  to 
565  nm.  Note  that  there  are  three  C2  Swan  bands  (d  ^flg  -  a  ^ITu  system)  present  in 


Table  4.7  Cl  Line  Intensities  at  725.66  nm  for  Solvents  Tested 


Solvent 

Intensity  (arb.  units) 

CCI4 

1.79x10^  ±603 

CHCI3 

8.69  X  10^^1440 

C2CI4 

3.32  X  10^^  +  306 

C2HCI3 

2.49  X  lO'^t  285 
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Intensity  (  X  1  0’  arb .  units)  Intensity  (  x  1  0^  arb .  units) 


600 


Figure  4.12  CCI4  Spectrum  Covering  585  -  685  nm  (100  pm  slit  width) 
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Figure  4.13  C2CI4  Spectrum  Covering  465  ■ 
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600 


-  565  nm 


the  spectrum.  (Spectra  for  CCI4  and  CHCI3  had  only  the  tiniest  hint  of  the  C2  Swan 
bands.)  They  arise  from  the  photodissociation  of  the  solvents  and  appear  as 
photofragment  emission  spectra.  Line  assignments  are  given  in  Table  4.8.'^^  Also 
notice  that  Figure  4. 14  depicts  a  C2HCI3  spectrum  very  similar  to  that  found  in 
Figure  4.13.  But  there  is  one  difference,  the  486. 1 3  nm  hydrogen  peak.^^’^^  At  first 
glance,  this  does  not  appear  to  be  seen  in  Figure  4. 1 3.  Nevertheless,  one  saw  the 
656.29  nm  hydrogen  peak  in  all  the  spectra.  It  would  make  sense  that  one  should 
also  see  the  486. 13  nm  peak  as  well.  In  fact,  it  is  there!  But  the  continuum 
background  is  so  large  that  it  masks  the  peak.  Remember  that  in  the  C2CI4 


Table  4.8  Spectrum  Line  Assignments  for  Figure  4.13 


Letter 

Element 

Wavelength 

v‘,v‘‘ 

A 

C2 

468.48  nm 

4,3 

B 

C2 

469.76  nm 

3,2 

C 

C2 

471.52  nm 

2,1 

D 

C2 

473.71  nm 

1,0 

E 

C2 

509.77  nm 

2,2 

F 

C2 

512.93  nm 

1,1 

G 

C2 

516.52  nm 

0,0 

H 

C2 

550.19  nm 

3,4 

I 

C2 

554.07  nm 

2,3 

J 

C2 

558.55  nm 

1,2 

K 

C2 

563.55  nm 

0,1 
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solution,  hydrogen  is  only  due  to  the  10  ppm  of  water  in  the  solution.  Since  the 
number  of  emitters  is  quite  low,  the  hydrogen  peak  most  likely  has  an  intensity 
much  lower  than  that  of  the  continuum  background. 

Figure  4. 15  depicts  a  CHCI3  spectrum  covering  the  wavelength  region  of 
345  to  445  nm.  Note  that  there  are  three  distinct  cyanogen  peaks  shown  (violet 
system,  B  Z  -  X  E).  Line  assignments  are  given  in  Table  4.9.  ’  Also  notice 
that  the  left  two  peaks  were  the  second  order  peaks  shown  in  Figure  4.11.  CCI4, 
C2CI4,  and  C2F1C13  had  similar  spectra  (see  Appendix  D). 


Figure  4.14.  C2HCI3  Spectrum  Covering  465  -  565  nm 
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Intensity  (  x  10^  arb.  units) 
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4.4  Saturated  Solvent  Solutions  (UV  Excitation) 

Saturated  solvent  solutions  were  run  at  a  10  }j,s  gate  width  and  a  1  )is  time 
delay.  In  every  case,  the  chlorine  bands  were  nonexistent.  Even  the  CN  bands 
were  missing.  This  can  be  seen  in  the  C2HCI3  spectrum  depicted  in  Figure  4. 16. 
(Similar  spectra  were  observed  for  CCI4,  CHCI3,  and  C2HCI3  as  well.)  One  will 
also  note  that  this  spectrum  is  veiy  reminiscent  of  the  deionized  water  spectrum 
shown  in  Figure  4. 1 7.  Line  assignments  are  noted  in  Table  4.10.  ’  The  one 
surprising  discovery  was  the  fact  that  nitrogen  appeared  in  every  case.  Now  could 
it  have  been  possible  that  the  gating  conditions  were  indeed  not  optimal?  As  a 
guide,  the  nitrogen  peaks  were  used  in  a  gating  experiment  (100  pm  slit  width). 
Spectra  were  gated  to  determine  an  optimal  delay  time.  Differing  gate  widths  were 
also  looked  at.  Sure  enough,  the  trend  showed  that  somewhere  around  a  1 0  ps 
gate  width  and  a  1  ps  time  delay  would  give  maximized  S/N  ratios.  Results  from 
the  gate  experiment  are  shown  in  Figure  4. 1 8  and  Table  4.11. 

On  a  final  note,  it  can  be  seen  from  Figures  4.16  and  4.17  that  intensities  in 
the  water  spectrum  are  about  four  times  greater  than  that  of  the  C2HCI3  spectrum. 
However,  the  water  spectrum  was  taken  with  a  100  pm  slit  width  (as  opposed  to 
50  pm).  This  will  take  care  of  a  factor  of  two.  So  when  slit  widths  are  taken  into 
account,  intensities  for  the  water  spectrum  double  that  of  the  C2HCI3  spectrum. 
Quenching  effects  (due  to  the  solvent)  may  be  one  possible  explanation  for  this. 
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Intensity  (  x  lO’arb.  units) 


Figure  4.16  C2HCI3  (1.20  g/l)  Spectrum  Covering  665  -  765  nm 


Figure  4.17  Water  Spectrum  Covering  665  -  765  nm  (100  pm  slit  width) 


46 


S/N  Ratio 


Table  4.10  Spectrum  Line  Assignments  for  Figures  4.16  and  4.17 


Table  4.11  Results  from  Gate  Width  Experiment  for  CCI4  (0.971  g/1) 


Gate  Width  (ns) 

Time  Delay  (ns) 

Signal 

Noise 

S/N  Ratio 

10 

0 

145 

10 

0.5 

330 

214 

10 

1 

202 

272 

10 

3 

41.8 

113 

10 

5 

88.0 

10.2 

8.63 

3 

1 

103 

267 

8 

1 

107 

208 

12 

1 

66.5 

314 

18 

1 

105 

253 
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V.  Discussion 


5.1  General 

As  can  be  seen  from  Figure  4.4,  the  continuum  background  due  to 
Bremsstrahlung  radiation  is  emitted  very  early  in  the  lifetime  of  the  plasma.  These 
results  indicate  (see  Table  4.2)  that  by  gating  the  signal,  one  can  significantly 
improve  the  S/N  ratio.  Note  that  the  analytical  lines  start  to  emerge  from  the 
continuum  background  and  successively  dominate  the  spectrum  as  the  time  delay 
is  increased.  Additionally,  Figures  4. 1  and  4.3  show,  that  at  later  times,  these 
analytical  lines  are  better  resolved.  This  is  due  to  the  fact  that  the  electron  density 
decreases  as  the  plasma  cools."*^ 

The  data  shown  in  Chapter  4  was  a  direct  result  of  surface  excitation  (as 
opposed  to  bulk  excitation)  of  the  liquid  sample  and  for  good  reason.  As  the  focus 
of  the  laser  waist  was  moved  further  into  the  bulk  of  the  solution,  the  spark  would 
start  to  fluctuate  and  eventually  disappear.  Consequently,  the  focus  was  kept  at  the 
surface  of  the  liquid  sample.  There  were  also  some  other  advantages  to  using 
surface  excitation.  First,  any  particles  suspended  in  the  solution  would  not  hamper 
spark  formation.  Second,  these  particles  would  not  scatter  light  to  any  great 
extent.^^ 
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One  might  recall  that  the  choice  was  made  to  run  the  laser  at  355  nm  shortly 
after  the  experiment  began.  This  was  done  to  look  at  the  wavelength  dependencies 
on  LIBS  signal  and  detection  limits.  First,  it  is  apparent  from  Figures  4.1  and  4.7 
that  the  LIBS  signal  is  much  greater  in  the  UV  than  in  the  IR.  This  made  sense  as 
absorption  processes  are  more  efficient  in  the  But  looking  at  Tables  4.2  and 
4.4,  it  can  be  seen  that  the  maximized  S/N  ratios  from  UV  and  IR  excitation  were 
virtually  the  same.  Even  the  detection  limits  were  very  close.  Initially,  this  seemed 
very  puzzling  as  the  LIBS  signal  in  the  UV  was  3  times  as  large  as  that  in  the  IR. 
But  looking  even  closer,  one  sees  that  the  noise  level  has  also  increased  in  the  UV. 
This  explains  the  almost  identical  S/N  ratios.  So  while  it  seems  that  UV  excitation 
will  give  a  larger  LIBS  signal,  there  seems  to  be  no  wavelength  dependence  on 
detection  limits. 

One  practical  reason,  however,  for  using  UV  excitation  would  be  to 
decrease  splashing  effects.  In  the  IR,  the  underside  of  the  prism  became  coated 
with  liquid  soon  after  the  start  of  the  experiment  thus  causing  the  laser  spark  to 
quickly  disintegrate.  However,  this  effect  was  much  less  pronounced  in  the 
ultraviolet.  The  reason  for  this  is  that  UV  excitation  is  more  efficient  than  IR 
excitation  in  the  sense  that  more  of  the  energy  goes  into  vaporizing  the  liquid  than 
heating  it.^ 
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5.2  Nickel  in  Solution 


In  looking  through  the  literature,  one  will  notice  that  there  are  very  few 
publications  available  on  the  laser-induced  breakdown  of  liquids.  However,  it 
might  be  instructive  to  compare  the  results  of  this  study  with  what’s  out  there. 

This  research  has  shown  detection  limits  for  nickel  in  solution  to  be  as  low  as 
29.4  mg/1  (in  the  UV).  Cespedes  et.  al.  also  looked  at  heavy  metals  in  solution  (Cd, 
Co,  Cr,  Cu,  Fe,  Mn,  Ni,  Pb,  Tl,  Zn  Au,  and  Hg).  Their  experimental  setup  was 
essentially  the  same  as  that  used  in  this  study,  but  for  most  of  the  metals  tested, 
their  results  yielded  far  worse  detection  limits.  In  fact,  a  detection  limit  of  1 500 
mg/1  was  obtained  for  Fe.  Some  metals,  such  as  Ni  and  Cr,  were  not  even 
detectable.  But  why  the  discrepancy?  The  most  likely  answer  for  this  is  the  fact 
that  Cespedes  et.  al.  used  bulk  excitation  of  the  solution  sample.^^  According  to 
Wachter  et.  al.,  surface  excitation  gives  far  better  detection  limits.  As  discussed  in 
Chapter  2,  their  study  lowered  detection  limits  for  uranium  in  solution  by  3  orders 
of  magnitude  (from  300  g/1  to  0. 1  g/1)  from  using  surface  excitation.  Detection 
limits  are  improved  by  surface  excitation  because  the  solution  is  vaporized  into  hot 
plasma  formed  above  the  liquid  surface  where  the  sample  is  more  efficiently 
excited.^^ 

Cremers  et.  al  looked  at  alkali  metals  in  solution  (Li,  Na,  K,  Rb,  and  Cs)  and 
found  detection  limits  generally  better  than  that  of  this  research  (0.006  -  1 .2  mg/1). 
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As  discussed  in  Chapter  2,  their  experimental  setup  differed  only  in  the  fact  that 
they  used  a  photomultiplier  tube  for  detection  purposes.^^  However,  this  was 
probably  not  the  reason  for  their  success.  The  fact  of  the  matter  is  that  alkali 
metals  are  easily  ionized  and  thus  should  provide  for  low  detection  limits.'^'^ 
Additionally,  it  should  be  pointed  out  that  2N/S  was  used  for  their  definition  of 
detection  limit.  This  alone  would  lower  their  results  significantly.^^ 

Knopp  et.  al  looked  at  LIBS  for  the  detection  of  metal  ions  (Cd^"^,  Pb^"^, 

Ba^"^,  Ca^^,  Li^"^,  Na"^,  Hg^^,  and  Er^^.  Since  ions  were  not  detected  in  this  study,  a 
direct  comparison  is  not  possible.  But  it  is  worth  taking  a  look  at.  As  mentioned 
in  Chapter  2,  their  experimental  setup  was  virtually  the  same  as  that  used  in  this 
research,  but  with  one  minor  exception.  Their  choice  of  excitation  was  an  excimer 
(A,  =  308  nm)  pumped  dye  laser  system.  Detection  limits  for  the  metal  ions  were 

generally  on  the  same  order  as  for  that  found  in  this  analysis  (except  for  Hg  and 

•2  1 

Er  which  were  not  detectable).  However,  as  with  Cremers  et.  al.,  their  lighter 
elements  (Ca^"^,  Li"^,  and  Na^)  had  much  lower  detection  limits  (in  the  pg/l  range).^^ 
Now  how  does  nickel  in  solution  compare  to  solid  nickel?  Thiem  et.  al. 
looked  at  solid  nickel.  Their  experimental  setup  differed  only  in  the  fact  that  they 
used  a  fiber  optic  bundle  to  collect  the  light  emitted  by  the  plasma.  They 
determined  a  detection  limit  of  10  ppm.  This  is  slightly  better  than  that  found  in 
this  study  (29.4  mg/1  ~  29.4  ppm).  Additionally,  Thiem  et.  al.  performed  their 
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experiment  in  an  ultra-high  vacuum.  If  anything,  this  helped  to  lower  their 
detection  limits  as  the  vacuum  would  reduce  the  background  signal  found  at 
atmospheric  pressure. 

On  a  final  note,  the  detection  limits  for  nickel  in  solution  were  found  to  be 
in  the  mg/1  range.  Hence,  this  technique  would  be  ideal  in  resolving  such  EPA 
compliance  issues  as  groundwater  contamination.  It’s  relative  simplicity  and  low 
cost  offers  significant  advantages  for  in  situ  monitoring  of  groundwater  runoff 
where  sample  preparation  may  not  be  an  option.  The  equipment  is  fairly  portable 
and  analysis  can  be  completed  within  15  to  20  minutes."*^ 

5.3  Pure  Solvent  Detection 

As  noted  in  Chapter  4,  the  Cl  line  intensities  at  725.66  nm  for  the  single 
bond  molecule  was  much  greater  than  that  of  their  double  bond  counterpart,  i.e., 
CCI4  >  C2CI4  or  CHCI3  >  C2HCI3.  Since  the  double  carbon  bond  is  stronger  than 
the  single  carbon  bond,  it  is  only  natural  that  it  would  take  more  energy  (heat)  to 
break  the  double  bond.  For  example,  as  seen  in  Appendix  E,  the  dissociation 
energy  for  C2CI4  is  approximately  3 1 9  kcal/mole,  while  for  CCI4,  it  is  3 1 1 
kcal/mole.  Hence,  it  is  not  surprising  that  the  Cl  line  intensity  for  the  single  bond 
molecule,  e.g.,  CCI4,  would  have  a  larger  peak  intensity  than  that  found  for  it’s 
double  bond  counterpart. 
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Also  mentioned  in  Chapter  4  was  the  fact  that  the  Cl  line  intensities  for  the 
chlorocarbons  were  much  higher  than  their  hydrocarbon  counterparts,  i.e.,  CCI4  > 
CHCI3  or  C2CI4  >  C2HCI3.  One  possibility  is  that  hydrogen  is  scavenging  some  of 
the  chlorine  atoms  to  form  hydrochloric  acid.  Unfortunately,  emission  lines  for 
HCl  were  not  detected.  Emission  for  HCl  usually  takes  place  in  the  ultraviolet 
(198-  237.5  nm)  where  it  is  often  hard  to  detect.  Pearse  and  Gayden  quote  that 
there  is  no  available  data  on  this  transition  but  it  is  likely  to  be  V  -  X 
Now  if  one  assumes  that  the  molar  excitation  efficiency  for  each  solvent  is  the 
same,  one  would  expect  that  the  chlorine  intensity  of  CHCI3  to  be  roughly  3/4  of 
CCI4.  If  this  assumption  is  true,  one  can  see  from  Table  4.7  (at  725.66  nm)  that 
CHCI3  should  have  an  intensity  of  approximately  1 .34  x  10^  ±  452.  However,  this 
estimate  is  too  high.  The  experimental  value  was  really  8.69  x  10'^  ±  440.  This 
leaves  a  difference  of  about  4.71  x  10'*.  This  discrepancy  is  most  likely  explained 
by  the  scavenging  of  Cl  atoms  spoken  of  earlier.  Likewise,  C2HCI3  should  have  an 
intensity  of  around  2.49  x  lO'^  ±  230.  Surprisingly,  it’s  experimental  value  was  this 
number.  Since  the  difference  is  zero,  one  notes  that  excitation  efficiencies  are 
probably  not  all  the  same  and  that  scavenging  of  the  Cl  atom  from  C2HCI3  is  really 
minute  at  best. 

There  is  one  simple  calculation  that  can  be  used  to  test  this  hypothesis. 

The  rate  law  for  the  amount  of  HCl  produced  is: 
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d  [HCl]  =  k  [Cl]  [H]  [M] 


(5.1) 


dt 


where  [M]  is  the  concentration  of  some  arbitrary  atom  used  to  advance  the  3 -body 
reaction.  This  equation  can  be  used  to  calculate  the  ratio  of  the  amount  of  HCl 
produced  from  CHCI3  to  that  from  C2HCI3.  Using  the  Cl  line  intensities 
(proportional  to  the  concentration  involved)  noted  above,  the  ratio  for  the  Cl  atoms 
is  8.69  X  IOV2.49  X  10"*  =  3.49.  Additionally,  the  ratio  for  the  H  atoms  (at  486.13 
nm)  is  4.72x10^1.35x10^  =3.50. 


Hence, 


d  [HCIJchcb  /dt  [Cl]cHCi3  [HJchcb 

d  [HC1]c2HC13  /dt  [C1]c2HC13[H]c2HC13 


«  (3.49)  (3.50) 
«  12.2 


This  means  that  CHCI3  is  approximately  12  times  more  likely  to  get  scavenged  for 
Cl  than  C2HCI3. 

There  is  one  other  way  to  determine  if  scavenging  is  going  on.  Let’s  start 
by  calculating  the  approximate  amount  of  chlorine  in  the  solvent.  The  following 
formula  may  be  used: 
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concentration  =  density  •  #  chlorine  ■  Mol.  Weight  ehlorine 


(5.2) 


per  solvent  Mol.  Weight  solvent 


The  results  of  this  equation  are  tabulated  in  Table  5. 1 .  Consider  CCU  and  CHCI3. 
One  notes  that  their  concentrations  are  within  10  %  of  each  other.  Assuming  their 
excitation  efficiencies  are  the  same,  one  would  expect  that  their  Cl  line  intensities 
should  be  within  10  %  of  each  other.  But  in  fact,  this  is  not  the  case.  One  can  use 
data  from  Tables  4.7  and  5. 1  to  generate  a  plot  of  intensity  versus  concentration. 
This  is  shown  in  Figure  5.1.  One  can  clearly  see  that  the  Cl  line  intensity  for  CCI4 
is  twice  that  for  CHCI3.  One  possible  explanation  for  this  is  the  scavenging  of  a  Cl 
atom  to  form  HCl. 

Now  one  can  also  use  equation  5. 1  to  look  at  the  concentration  of  carbon  in 
the  solvent.  Substituting  for  earbon  into  equation  5.1,  one  gets  the  results  listed  in 
Table  5.2.  As  expected,  one  notes  that  the  double  carbon  bond  has  a  carbon 
concentration  of  approximately  twice  that  of  its  single  bond  counterpart.  Now  one 
can  also  use  the  data  from  Tables  5.2  and  5.3  to  generate  a  plot  of  intensity  versus 
eoncentration.  (CN  intensities  can  be  used  because  carbon  can  only  come  from  the 
solvent  and  for  every  carbon  atom  emitted,  there  is  one  nitrogen  atom  emitted.) 
This  plot  is  shown  in  Figure  5.2.  Considering  the  fact  that  CCI4  and  CHCI3  have 
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concentrations  within  approximately  1 5  %  of  each  other,  one  would  expect  that 
their  intensities  would  be  very  close  to  each  other.  And  this  is  the  case.  But  when 
one  looks  at  C2CI4  and  C2HCI3,  there  is  a  different  story  to  tell.  C2HCI3  has  a  much 
higher  intensity  than  that  of  C2CI4.  The  author  is  not  sure  why  this  is  happening. 
Further  investigation  will  hopefully  answer  this  mystery.  Nevertheless,  using  the 
data  from  Tables  4.7  and  5.3,  one  can  calculate  Cl/CN  intensity  ratios.  This  is 
shown  in  Table  5.4.  One  notes  that  the  ratios  for  the  chlorocarbons  are 
approximately  twice  that  of  the  hydrocarbons.  This  would  also  suggest  that 

o  1 

scavenging  of  Cl  is  going  on. 


Table  5.1  Concentration  of  Chlorine  in  the  Solvent  Tested 


Solvent 

Density  of 
Solvent  (g/ml)'^* 

Mol.  Weight  of 
Solvent"** 

Cl  Concentration 
(K/1) 

CCI4 

1.58 

154 

CHCI3 

1.48 

119 

C2CI4 

1.61 

166 

1.46 

131 

■nEHQHi 
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Cl  Intensity  (  x  lO’arb.  units) 


Figure  5.1  Cl  Intensity  Versus  Cl  Concentration  at  725.66  nm 


Table  5.2  Concentration  of  Carbon  in  the  Solvent  Tested 


Solvent 

Density  of 
Solvent  (g/ml)'^* 

Mol.  Weight  of 
Solvent'*^ 

C  Concentration 

(g/l) 

ecu 

1.58 

154 

123 

CHCU 

1.48 

119 

149 

C2CI4 

1.61 

166 

233 

C2HCI3 

1.46 

131 

267 
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CN  Intensity  (x  1  0^  arb.  units) 


Table  5.3  CN  Line  Intensities  at  359.04  nm  for  Solvents  Tested 


120  160  200  240  280 


C  Concentration  (g/1) 

Figure  5.2  CN  Intensity  Versus  C  Concentration 
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Table  5.4  CI/CN  Ratios  for  Solvents  Tested 


Solvent 

Cl/CN  Ratio 

ecu 

10.2 

CHCls 

5.23 

C2CI4 

2.68 

C2HCI3 

1.06 

Let’s  revisit  equation  5.1.  Bauch  et.  al.,  using  argon  as  the  third  body, 
expressed  the  rate  constant  as: 

k  =  2.60  X  1 0*^  exp  (9940  /  T)  cm^  mole’^  s'*  (5.3) 

Assuming  a  temperature  (upon  recombination)  of  1,000  K,  the  rate  constant 
k  =  5.39  X  10*^  cm^  mole'^  s'* «  1.49  x  10'^**  cm^  molecule'^  s'*.  Additionally, 
equation  5. 1  can  be  written  as  a  pseudo-second  order  equation.  If  the 
concentration  of  argon  is  kept  constant, 

d  [HCl]  «  k’  [H]  [Cl]  (5.4) 

dt 
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where  k’  =  k  [Ar],  Using  the  ideal  gas  law  (p  =  NkT)  with  p  =  1  atm  and  T  = 
1,000  K,  N  «  [Ar]  =  7.34  x  10‘*  cm'^  Hence,  k’  «  1 .09  x  10'"  cm^  molecule’^  s’^ 
As  one  can  see,  this  is  a  fairly  fast  rate  constant  and  recombination  should  occur 
rather  quickly.  This  lends  credence  to  the  idea  that  H  is  scavenging  Cl  to  form 

j^Ci3,31,46 

Back  in  Chapter  4,  it  was  revealed  that  CCI4  and  CHCI3  had  only  the  tiniest 
hint  of  the  C2  Swan  bands.  This  should  come  as  no  surprise.  When  C2CI4  and 
C2HCI3  dissociate,  there  will  be  excess  C2  molecules  around  to  form  the  C2  Swan 
bands.  However,  when  CCI4  and  CHCI3  dissociate,  there  are  only  C  atoms 
available,  not  C2.  To  form  the  C2  Swan  bands  from  CCI4  and  CHCI3,  these  C 
atoms  must  recombine  to  form  C2.  The  chance  of  this  happening  is  relatively  low 
because  carbon  is  most  likely  to  recombine  with  H  (Cl)  to  form  CH  (CCl),  i.e.,  the 
heats  of  formation  (enthalpy)  are  A  Hcci  =111  kcal/mole,  A  Hch  =142  kcal/mole, 
and  A  Hc2  =  200  kcal/mole."^^  Unfortunately,  neither  CH  or  CCl  were  detected  by 
this  technique. 
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5.4  Saturated  Solvent  Detection 


As  seen  in  Chapter  4,  saturated  solvent  solutions  were  not  detectable  with 
LIBS.  They  basically  mirrored  the  water  spectrum.  However,  there  is  a  way  to 
determine  if  saturated  solvent  detection  was  even  possible  in  the  first  place.  Using 
data  from  Table  4.7  along  with  the  origin  as  a  data  point,  one  can  generate  a 
calibration  curve  for  each  solvent  in  question.  This  is  shown  in  Figures  5.3 
through  5.6.  Detection  limits  can  be  calculated  using  equation  4. 1  and  are  shown 
in  Table  5.5.  Looking  back  in  Ch.  3,  one  can  see  that  the  detection  limits  for  the 
solvents  (chlorine)  in  question  are  much  higher  than  their  solubility  limits.  Hence, 

Q 1 

one  would  not  be  able  to  detect  the  solvent  (chlorine)  in  the  saturated  solution. 


Figure  5.3  Simulated  Calibration  Curve  for  CCI4 
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Concentration  (g/1) 


Figure  5.5  Simulated  Calibration  Curve  for  C2CI4 


Figure  5.6  Simulated  Calibration  Curve  for  C2HCI3 


Table  5,5  Calibration  Curve  Data  for  Chlorine 


Solvent 

Noise  N 

CCI4 

WMBaaSMM 

123 

39.5 

CHCI3 

HBEIEDS9Hi 

65.8 

46.0 

C2CI4 

588 

24.1 

73.2 

■■ssshb 

588 

20.9 

84.4 

64 


VI.  Conclusions  and  Recommendations 


This  study  has  presented  the  spectroscopic  technique  known  as  Laser- 
Induced  Breakdown  Spectroscopy.  It  has  been  applied  to  nickel  in  solution.  Limits 
of  detection  ranged  between  29.4  and  56.1  mg/1.  Detection  limits  of  the  sort 
should  provide  a  moderate  degree  of  analytical  performance  for  such  EPA 
compliance  issues  as  groundwater  monitoring.  LIBS  was  also  applied  to  solvents 
such  as  CCI4,  CHCI3,  C2CI4,  and  C2HCI3.  Cl,  C2,  and  CN  bands  were  detected  in 
the  pure  solvent  solutions.  However,  solvent  detection  was  not  possible  for  the 
saturated  solutions.  Continued  work  with  pure  solvent  solutions  may  lead  to  future 
work  with  saturated  solvent  solutions. 

As  with  any  experiment,  there  probably  are  some  improvements  that  can  be 
made  to  the  equipment.  The  following  recommendations  are  suggestions  that 
might  even  help  lower  detection  limits. 

1 .  A  Nd;  YAG  laser  has  served  as  the  main  source  of  excitation  throughout 
the  entire  course  of  this  experiment.  While  it  has  fulfilled  it’s  role,  the 
addition  of  a  second  laser  could  be  used  to  increase  emission  from  the 
plasma  and  thereby  lower  detection  limits.  Once  the  plasma  has  been 
initiated  by  the  first  pulse,  a  second  pulse  would  be  used  to  further  excite 
the  plasma.  However,  the  second  pulse  would  only  interact  with  the  plasma 
above  the  sample.  Of  course,  necessary  steps  would  have  to  be  taken  to 
ensure  synchronization  of  the  two-laser  system.^^ 

2.  Optical  alignment  of  the  laser  axis  is  the  key  to  good  quantitative  results. 
The  use  of  focusing  lenses  and  2-D  translation  stages  to  perform  this  task 
was  difficult  at  best.  One  might  instead  use  a  fiber  optic  cable  both  to 
deliver  laser  radiation  to  the  solution  sample  and  to  collect  light  emitted 
from  the  focal  volume. 
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3.  In  this  study,  detection  of  the  spark  plasma  was  triggered  off  of  the 
Nd.YAG  laser.  Unfortunately,  the  firing  of  the  laser  did  not  always  result 
in  breakdown  of  the  liquid  Therefore,  an  optical  trigger  system  may  be 
more  appropriate  for  this  application.  Signal  processing  would  only  take 
place  when  the  laser  spark  was  observed.  An  optical  trigger  detector  is 
available  commercially  through  EG&G  PARC."^ 

4.  For  those  who  are  not  concerned  with  in  situ  detection,  one  might  try 
performing  this  experiment  in  a  vacuum.  This  should  reduce  the 
background  signal  found  at  atmospheric  pressure  and  thus  lower  detection 
limits. 

5.  Saturated  solvent  detection  proved  to  be  uneventful  in  this  experiment. 
However,  one  might  try  running  the  experiment  in  reverse,  starting  with  the 
pure  solvent.  Next,  a  small  amount  of  water  would  be  added  to  the  solvent 
and  the  experiment  would  be  run  again.  This  would  continue  until  the 
solvent  itself  was  saturated.  Any  change  in  spectral  features  (from  one 
solution  to  the  next)  could  be  indicative  of  the  amount  of  solvent  in  the 
solution.^ 


Appendix  A  -  Area/Noise  Ratios 


The  area/noise  ratio  is  readily  calculated.  First  of  all,  the  area  is  simply  the 
area  under  the  curve  (peak).  A  nice  software  package  can  easily  perform  the 
integration.  For  this  study,  the  361 .94  nm  peak  for  Ni  (UV)  was  used.  Next,  noise 
is  obtained  from  a  portion  of  the  Ni  spectrum  that’s  essentially  featureless.  For  this 
region,  the  intensities  and  wavelengths  are  known  so  that  the  multiplication  of  the 
two  give  rise  to  an  area  for  the  noise.  Finally,  take  the  standard  deviation  of  the 
noise  values.  Once  completed,  one  can  easily  calculate  the  area/noise  ratios  as 
seen  in  the  table  below.  Units  for  area  and  noise  are  in  arb.  units  -  nm.'*^ 


Table  A.l  Area/Noise  Ratios  for  Ni  (UV)  at  361.94  nm 


Gate  Width  (^is) 

Time  Delay  (jis) 

Area 

Noise 

A/N  Ratio 

10 

1 

15.1 

10 

_ 

3 

IliEBEBSliH 

18.8 

10 

6 

HBESBSHl 

14.9 

10 

9 

3.62 

10 

12 

■BOSH 

111 

3 

0.5 

miE^ssysi^^i 

4.04 

7 

0.5 

5.58 

10 

0.5 

ggg^ggi 

10.3 

14 

0.5 

8.91 

18 

0.5 

HIKIHill 

gg^^gigi 

7.39 
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Appendix  B  -  Systematic  Error 


One  source  of  error  in  the  experiment  is  systematic  error.  They  cause  the 
results  of  replicate  measurements  to  be  either  high  or  low.  There  are  three  types  of 
systematic  error.  The  first  type  is  instrument  error.  They  are  caused  by 
imperfections  in  the  measuring  devices  used  in  the  experiment.  For  example, 
pipets  and  volumetric  flasks  have  volumes  slightly  different  from  that  indicated  by 
their  graduations.  The  next  type  of  error  is  the  method  error.  They  arise  from  the 
nonideal  chemical  behavior  of  the  reagents  upon  which  analysis  is  based.  Sources 
of  nonideality  include  instability  of  the  species  tested  and  the  possible  occurrence 
of  side  reactions  that  interfere  with  the  measurement  process.  Finally,  the  last  type 
of  systematic  error  is  personal  error.  They  result  from  the  carelessness  and 
inattention  of  the  experimenter.  Additionally,  personal  error  can  arise  from  bias  on 
the  part  of  the  observer.  Most  people  have  a  natural  tendency  to  estimate 
measurements  in  such  a  way  as  to  improve  their  results.  The  author  of  this  study 
suspects  that  the  systematic  error  in  this  experiment  probably  outweighs  the 
statistical  error  inherent  in  the  results. 
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Appendix  C  -  Determination  of  Noise  (N)  at  Various  Concentrations 


Noise  is  calculated  by  the  following  means.  For  a  portion  of  the  Ni 
spectrum  that  is  essentially  featureless,  intensity  values  are  known.  N  is  the 
standard  deviation  of  those  values.^  Units  for  N  are  in  arbitrary  units. 


Table  C.l  Noise  Levels  at  Various  Concentrations 


Gate  Width  (ps) 

Time  Delay  (jis) 

Excitation 

Concentration 

(mg/1) 

N 

50 

8 

IR 

1000 

8.38 

50 

8 

IR 

200 

12.0 

50 

8 

IR 

100 

10.4 

10 

3 

UV 

1000 

20.9 

10 

3 

uv 

500 

17.2 

10 

3 

UV 

200 

27.4 

10 

3 

uv 

125 

14.6 

10 

3 

uv 

100 

11.7 

10 

3 

uv 

62.5 

18.9 

10 

3 

uv 

50 

20.0 
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Appendix  D  -  Spectra  for  Solvents  Tested 


Displayed  in  this  appendix  is  a  series  of  40  spectra  taken  of  the  four 
solvents  studied  in  this  research.  All  spectra  are  taken  at  a  50  pm  slit  width  unless 
otherwise  noted.  Figures  D.  1  through  D.8  show  two  chlorine  peaks  at  725.66  and 
754.71  nm.  These  figures  also  depict  second  order  cyanogen  bands  as  well  as 
nitrogen  lines.  Line  assignments  are  listed  in  Table  4.6.  Figures  D.9  through  D.16 
show  the  656.29  hydrogen-alpha  peak.  Figures  D.17  through  D.28  introduce  the 
C2  Swan  bands.  Line  assignments  are  given  in  Table  4.8.  Note  that  Figures  D.26 
and  D.28  picture  the  486. 13  nm  hydrogen  peak.  Finally,  Figures  D.29  through 
D.40  show  the  cyanogen  bands.  Line  assignments  are  listed  in  Table  4.9. 
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Intensity  (  x  1  O'  arb.  units)  Intensity  (  x  10'  arb.  units) 


Figure  D.6 


ngth  (nm  ) 


ng  665  -  775  nm  (100  fim  slit  width) 


Intensity  (  x  1  0^  arb.  units) 


Figure  D.7  C2CI4  Spectrum  Covering  665  -  775  nm 


Figure  D.8  C2HCI3  Spectrum  Covering  665  -  775  nm 


74 


Intensity  (  x  1  0^  arb.  units) 


Figure  D. 


Figure  D.IO 


Intensity  (  x  1  0^  arb.  units)  Intensity  (  x  1  0^  arb.  units) 


Figure  D.ll  C2CI4  Spectrum  Covering  625  -  725  nm 


Figure  D.12  C2HCI3  Spectrum  Covering  625  -  725  nm 
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Intensity  (  x  10^  arb.  units) 


Figure  D.14 


Intensity  (  x  lO’arb,  units)  Intensity  (  x  lO^arb.  units) 
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Figure  D.16  C2HCI3  Spectrum  Covering  585  -  685  nm 
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Intensity  (  x  1  0^  arb  units)  Intensity  (  x  1  0’  arb.  units) 


Figure  D.18 


Intensity  (  x  10’  arb.  units) 


Figure  D.19  C2CI4  Spectrum  Covering  545  -  645  nm 


Figure  D.20  C2HCI3  Spectrum  Covering  545  -  645  nm 


0 


8i 


Intensity  (  x  10^  arb.  units) 


Figure  D.21  CCI4  Spectrum  Coveriog  505  -  605  nm  (100  pm  slit  wid 


Figure  D.22  CHCI3  Spectrum  Covering  505  -  605  nm  (100  pm  slit  wu 
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Intensity  (  x  1  0^  arb  units)  Intensity  (  x  1  0^  arb.  units) 


400 


Wavelength  (nm) 


Figure  D.27  C2CI4  Spectrum  Covering  465-565  nm 


Figure  D.28  C2HCI3  Spectrum  Covering  465  -  565  nm 


8. 


Intensity  (  x  1  0^  arb,  units)  Intensity  (  x  1  O'  arb.  units) 


Figure  D.39  C2CI4  Spectrum  Covering  345  -  445  nm  (0.8  ND  Filter) 


Figure  D.40  C2HCI3  Spectrum  Covering  345  -  445  nm  (0.8  ND  Filter) 
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Appendix  £  -  Bond  Dissociation  Energy 


The  bond  dissociation  energy  for  a  molecule  can  be  calculated  from  the 

o  1 

heats  of  formation  (enthalpy)  AH  according  to  the  following  formula  : 


^^Hreaction  ^  ^  '^Hj-eagents 


(E.l) 


For  the  reaction  (using  data  from  Table  E.l):  CCI4  ->  C  +  4C1 


AHreaction  =  311  kcal/molc 


For  the  reaction  (using  data  from  Table  E.  1 );  C2CI4  C2  +  4C1 


AHreaction  =  319  kcal/molc 


Table  E.  1  Heats  of  Formation  for  Several  Atoms/Molecules 


Atom/Molecule 

Heat  of  Formation 
(kcal/mole)'*^ 

C 

172 

Cl 

29.0 

C2 

200 

CC14 

-22.9 

C2C14 

-2.70 
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